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Abstract. We investigate the moments of the derivative, on the unit circle, of charac- 
teristic polynomials of random unitary matrices and use this to formulate a conjecture for 
the moments of the derivative of the Riemann £ function on the critical line. We do the 
same for the analogue of Hardy's Z-function, the characteristic polynomial multiplied by a 
f-H ' suitable factor to make it real on the unit circle. Our formulae are expressed in terms of a 

determinant of a matrix whose entries involve the I-Bessel function and, alternately, by a 
combinatorial sum. 

!: 

1. Introduction 

Characteristic polynomials of unitary matrices serve as extremely useful models for the 
Riemann zeta-function ((s). The distribution of their eigenvalues give insight into the dis- 
' tribution of zeros of the Riemann zeta-function and the values of these characteristic poly- 
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nomials give a model for the value distribution of C( s )- See the works [KS] and [CFKRS] for 
detailed descriptions of how these models work. The important fact is that formulas for the 
moments of the Riemann zeta-function are suggested by the moments of the characteristic 
^ ■ polynomials of unitary matrices. 

We consider two problems here: the moments of the derivative of the characteristic poly- 
nomial Aa(s) of an N x iV unitary matrix A, and also the moments of the analogue of 
the derivative of Hardy's Z-function, the characteristic polynomial multiplied by a suitable 
factor to make it real on the unit circle. 

In its simplest form our problem is to give an exact formula, valid for complex r with 



> 



X 



9ftr > 0, of the moments of the absolute value of the derivative of characteristic polynomials 

(1.1) / \A' A (l)\ r dA N 

JU(N) 

or of 

(1.2) / \Z' A {\)\ r dA N . 

JU(N) 

Here we are integrating against Haar measure on the unitary group, and Za{s) is equal to 
Aa(s) times a rotation factor that makes it real on the unit circle. See the next section for 
the precise definition. 
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Unfortunately, we cannot yet solve either of these problems. However, we can give as- 
ymptotic formulas when r = 2k for positive integer values of k. The first two of these 
formulas involve the Maclaurin series coefficients of a certain kxk determinant, while the 
third involves a combinatorial sum. 

Theorem 1. For fixed k and N — > oo we have 



(1.3) / \A' A (l)\ 2k dA N ~b k N k ' +2k , 

Ju(N) 

where 

(1.4) K = (-l)^t (*) (ff + \e-'z-^ 



x=0 



and I v (z) denotes the modified Bessel function of the first kind. 
Theorem 2. For fixed k and N — > oo we have 

(1.5) / |^(l)| 2fc ^~6^iV fc2+2fc , 

-/[/(AO 



where 

2k 



1.6) *4 = (-l)^V 2 ( A) (^e-f x" fc2 / 2 det (/^■_ 1 (2V5)) 



x=0 



We also have combinatorial description of b' k . 
Theorem 3. 

(1-7) 6t = 

(-i) fe(fe+1)/2 E me p- 1(2fe) © r (nL (ilw* j)) , 

where Pp +1 (2k) denotes the set of partitions m = (m , . . . , m k ) of 2k into k + 1 non-negative 
parts. 

We have computed some values of b k and b' k ; these are tabulated at the end of the paper. 
Applying the random matrix theory philosophy suggests the conjecture: 

Conjecture 1. 

(1.8) \ jf |C'(l/2 + it)\ 2k dt ~ a k b k log(T) fc2+2fe 
and, similarly for Hardy's Z function, 

(1.9) i jT |Z'(l/2 + ft)! 2 ** ~ a k b' k log(T) fe2+2fe , 
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where a k is the arithmetic factor 
(1.10) a k = Tl 



p) m\Y(k) I P 

m=0 v v ' 



Remarks: 

(1) The factor a k is the same arithmetic contribution that arises in the moments of the 
Riemann zeta function itself, see [KS] or [CFKRS]. For an explanation of why these 
moments factor, asymptotically, into the product of a contribution from the primes, 
dfc, and a coefficient calculated via random matrix theory, see [GHK]. 

(2) In this paper we are only concerned with the leading asymptotics of the moments of 
C'(l/2 + it) and Z' (1/2 + it). Consequently, we use the fc-fold integrals for moments 
given below in Lemma El If one wishes to study lower order terms one would need to 
use the full moment conjecture for ( and Z given in [CFKRS] as a 2k fold integral. 

(3) Forrester and Witte have taken our Theorems 1-2 and managed to find an alternate 
expression for b k and b' k involving a Painleve III' equation, and also an expression 
involving a certain generalised hypergeometric function [FW, section 5]. 

(4) In his PhD thesis, Chris Hughes gives a similar conjecture for a more general problem 
involving mixed moments [Hug, conj 6.1]. What is new in our paper are the formulas 
for bk and b' k . For comparison, we state Hughes' formulation of the conjecture for 
Hardy's Z function. Let 

I{h,k) = / \Z{t)\ 2k ~ 2h \Z'{t)\ 2h alt. 



Hughes conjectures that 

I(h,k) ~ B(h,k)a k f k T(logTf +2h 
where a k is given above, 



fc-1 .. 



and B(h, k) is the constant that Hughes obtains from the analogous moments for Za 
via random matrix theory: 

2h 



" n=0 ^ ' 

where 

(2k — n + i — 1)! f i + k - n - I + m\ f i + m - l\ 
l ' j ~ ^ (2k-n + i-l + m)\ \ m )\ j -1 J 

To get Hughes' conjecture in this form, see (6.60), (6.51), and (6.52) of [Hug] and 
replace j3 by (3/{iN). 
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(5) Brezin and Hikami [BH] attempt to use a similar approach to obtain a theorem for 
moments of derivatives of characteristic polynomials, but there is an error in their 
paper. 

(6) Numerically, we have observed that bk ~ 4 _ as k — > oo. Other than the power of 
4, the r.h.s is the constant that appears in the moments of characteristic polynomials 
[KS]: 

/ \A A {l)\ 2k dA N ~ f k N k \ asiV^oo. 

Ju(N) 

A heuristic explanation is as follows. The large values of |A^(1)| occur near the large 
values of (A^l)!, namely when all the eigenvalues are close to —1. The derivative 
of Aa involves a sum of N terms, each of which is missing, when the eigenvalues are 
close to —1, one factor of size roughly 2. A comparison with the 2/c-th moment of 
jA^I thus gives an extra N 2k and a factor of 2~ 2k . We have not attempted to make 
this arguement rigorous. We have also not attempted to show that b k and b' k are 
non-zero. 

(7) The problem of moments of the derivative, are related, through Jensen's formula, 
to the problem of zeros of the derivative. This approach requires knowledge of the 
complex moments of the derivative and we are only able to obtain integer moments. 
For characteristic polynomials one is interested in studying the radial distribution of 
the zeros of the derivative. Francesco Mezzadri has the best results in this direction 
[Mez]. On the number theory side, one is interested in the horizontal distribution 
of the zeros of Partial results have been obtained by Levinson and Montgomery 
[LM], Conrey and Ghosh [CG], Soundararajan [Sou], and Zhang [Z]. 



2. Notation 

If A is an A" x A" matrix with complex entries A = we let A* be its conjugate 

transpose, i.e. A* = (bjk) where b^ = OkJ. A is said to be unitary if AA* = I. We let U(N) 
denote the group of all A^ x A" unitary matrices. This is a compact Lie group and has a Haar 
measure. 

All of the eigenvalues of A e U(N) have absolute value 1; we write them as 
(2.1) e ie \e ie \...e ldN . 

The eigenvalues of A* are e~ t01 , . . . , e~ t0N . Clearly, the determinant, det A = Yln=i e%dn °f a 
unitary matrix is a complex number with absolute value equal to 1. 

We are interested in computing various statistics about these eigenvalues. Consequently, 
we identify all matrices in U(N) that have the same set of eigenvalues. The collection of 
matrices with the same set of eigenvalues constitutes a conjugacy class in U(N). Weyl's 
integration formula [Weyl, pg 197] gives a simple way to perform averages over U(N) for 
functions / that are constant on conjugacy classes. Weyl's formula asserts that for such an 
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(2.2) / f(A) cfflaar = / f(6 1 , 6 N )dA N 

JU(N) J[0,2n] N 

where 

(2.3) dA N = J] | e * ^d6 l ...d6 N 



JV!(2tt) 



e J i 

i<j<k<N 

The characteristic polynomial of a matrix A is denoted A A (s) and is defined by 

AT 

(2.4) A A (s) = det(J - sA*) = JJ(1 - se~ i6n ). 

n=l 

The roots of A A (s) are the eigenvalues of A and are on the unit circle. Notice that this 
definition of the characteristic polynomial differs slightly from the usual definition in that it 
has an extra factor of det(A*). We regard A A (s) as an analogue of the Riemann zeta-function 
and this definition is chosen so as to resemble the Hadamard product of (. 
The characteristic polynomial satisfies the functional equation 

N N 

(2.5) A A (s) = (- a )"J] e -*»J](i- e *»/ a ) 

n=l n=l 

(2.6) = (-lfdetA* s N A A .(l/s). 
We define the ^-function by 

(2.7) Z A (s) = e - mN / 2 e^^ 2 s- N / 2 A A ( S y, 

here if N is odd, we use the branch of the square-root function that is positive for positive 
real s. The functional equation for Z is 

(2.8) Z A (s) = (-l) N Z A .(l/s). 
Note that 



(2.9) Z A (e») = Z A (e* & ) 

so that Z A (e td ) is real when 9 is real. We regard Z A (e ld ) as an analogue of Hardy's function 
Z(t). 

We let I n be the usual modified Bessel function with power series expansion 



(2.io) ux) = (-) n y 

The way that the I-Bessel function enters our calculation is through the following formula: 
2niJ ]zl=1 z 2 * dZ ~ (Lt)^/ 2 ■ 
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This formula can be proven by comparing the coefficient of z 2k ~ 1 in e Lz+t ^ z with the power 
series formula for l<ik-\- 

We let A(zi, . . . , z k ) denote the Vandermonde determinant 

(2.12) A(z 1 ,...,z k ) = det(zr 1 ). 

kxk 

We often omit the subscripts and write A(z) in place of A(z±, . . . , z k ). Also, we allow 
differential operators as the arguments, such as 

The key fact about the Vandermonde is that 

(2.14) A{z u ...,z k ) = 11 (Zj-Zi). 

We let 

(2.15) z(x) = -±— = 1 + 0(1). 

1 — e x x 

The function z(x) plays the role for random matrix theory that ((1 + x) plays in the theory 
of moments of the Riemann zeta-function. See for example pages 371-372 of [CFKRS2]. 
We let S denote the subset of permutations a e of {1,2,..., 2k} for which 

(2.16) cr(l) < <t(2) < • • • < a{k) 
and 

(2.17) a(k + 1) < a(k + 2) < ■■■ < a {2k). 

We let Pq +1 (2/c) be the set of partitions m = (m , . . . , m k ) of 2k into k + 1 non-negative 
parts. This quantity arises from the multinomial expansion 

(2.18) (so + *i + • • • + a*) 2 * = E ( 2 £) x °°"- xmi 



where 
(2.19) 



mdP^ +l (2k) 

2k\ [2k)\ 



k 



m J mo'. . . . m k \ 
3. Lemmas 



The main tool in proving theorems 1-2 is to take formulas (Lemma 3) for moments of 
characteristic polynomials with shifts, differentiate these with respect to the shifts, and then 
set the shifts equal to zero. This gives fc-fold contour integrals. To separate the integrals 
involved, we introduce extra parameters and differential operators to pull out a portion of 
these integrands. 
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Lemma 1. Assume that aci, . . . , ct2k are distinct complex numbers. We have 

(3.1) f JjA A (e-°«)A x .(e a ^*) dA N = ^ e ^£*=iK(;)-«;) JJ z (a a{j) - a a{k+i) ). 

Ju ( N ) j=l aeS l<i,j<k 

This is proven in section 2 of [CFKRS2]. See formulas (2.5), (2.16), and (2.21) of that 
paper. The definition given there of the characteristic polynomial differs slightly from the 
one we use here, and that introduces some extra exponential factors in (2.21) of the afore- 
mentioned paper, and also necessitates replacing the a's by —a's. 

Since 

(3.2) Z A (e- a i)Z A *(e a >+ k ) = {-1) N e N ^- a ^ k)/2 A A {e~^)A A *{e^+ k ) 
we can write a corresponding lemma for Z. 

Lemma 2. Assume that ati, . . . , a2k are distinct complex numbers. Then 

k 



(3.3) / \[Z A {e-^)Z A ,{e^ k ) dA t 

JU(N) =1 



N . 



We can express the sums in the last two lemmas as integrals. Thus we have 
Lemma 3. Assume that all of the ctj are smaller than 1 in absolute value. Then 



(3.4) / T\A A (e- a i)A A *(e a >+ k ) dA N 

JU(N) fJl 

I e ^E*UK-%0 Y[ ziwi-a^Wziwi-w^Wdw 

J \w;\ = l i^.-^i- s_l: _• 1 



k\(2ni 



i<i<k i+j j=\ 

l<j<2k rJ 



and 



(3.5) / f[Z A (e-^)Z A *(e a ^) dA N 

JU(N) j=1 

— — J] 2fe a k 

= (- i )" fc£ ik^ / ^ Uwi n ^ - «i) - ^m^r 

J\ Wi \=i , <i<k ... ._, 



i<i<k i+j j=\ 

l<j<2k ^ J 



In this lemma, and its corollary below, we do not require the a's to be distinct. The proof 
of Lemma 01 is a straight-forward evaluation of the residues in the integral in ()3.4|) . arising 
from the factor z(wi — a,), to obtain the ( 2 fc fc ) terms in (|3. lj) . Each of the k integrals in ()3.4|) 
results in a sum over 2k residues, but due to the factor Yli^j z ( w i ~ w j)~ 1 ^ an y one °f these 
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2k 2 terms is zero if the residue of two of the integrals, say Wi and Wj, are evaluated at the 
same point ot£. 

Using the fact that z(w) = 1/w + 0(1) we easily deduce 

Corollary 1. Suppose that otj = (Xj(N) and \otj\ <C 1/N for each j . Then 
(3.6) / TTA A (e-^)AA*(e a ^) dA N 

JU(N) j=1 

= \ / e^E^iK-aj) n ¥j -(wi-w 3 -) TJdw +Q(N k2 - 1 ) 



with an implicit constant independent of N; similarly, 
k 



(3.7) / T\Z A (e-^)Z A ,(e a ^) dA r 

JU(N) j=1 



— — \^ 2fc rv TT / \ k 



Lemma 4. Let f be k — 1 ^mes differentiable, k > 1. TTien 

^)fl/W=det (,«->(!,)) 

where by A(d/dL) we mean the differential operator 
. . T-r / d d \ , / d^ 1 



v dL,- dLj / kxk\d]j. 1 

Proof. This follows using the definition of the Vandermonde determinant. Noticing that 
row i of the matrix only involves L iy we factor the product into the determinant. 

Lemma 5. Let f be 2k — 2 times differentiable. Then 

(^) A 2 (^)(lW.>)| = *!det(/™"<>(L)). 

More generally, suppose that g(L 1 , . . . , L&) = X]f=i a r Ili=i fr,i(Li) is a symmetric function 
of its k variables. Then 

(3.10) A 2 (A) g{Ll) . . . , Lk ) |^ ^ = fc! £ det (i^ 2) (L)) . 



MOMENTS OF THE DERIVATIVE 9 

Proof. Applying the Vandermonde a second time to Lemma 0] we get 

< 3 - u > K£>- (/U " ,(i,)) ' 

Expand the determinant as a sum over all permutations \i of the numbers 1,2, ... ,k: 

k 

(3.12) det(/^(L0) =X>gn(/x)n/ W-1 (^)- 

fi i=l 

Apply Lemma 0] to find that a typical term above equals 

(3.13) sgn(/i)det(/^- 2 )(L i )). 

kxk 

Setting Li = L for 1 < i < k, we may rearrange the rows so as to undo the permutation \x. 
This introduces another sgn(/x) in front of the determinant and gives 

(3.14) det(f i+j - 2 \L)). 

kxk 

Since there are k\ permutations /i, we get 

(3.15) k\det(f {i+j - 2) (L)). 

The proof of the second part of the lemma is left to the reader. 

Lemma 6. Suppose that P andQ are polynomials with Q{w) = Ylj=i( w ~ a j) an ^ niax \atj\ < 
c. Then 



Proof. Since 



(3.17) p(4t )e wL = e wL P(w) 



\alL / 

the derivatives can be pulled outside the integral immediately. With the Laplace transform 
pair e xa and — , related by 

(3.18) e xa = / dw, 

2-ni J\ w \ =c w — a 

we merely apply repeatedly the Laplace convolution formula, which for Laplace transform 
pairs fi and fa states that 

(3.19) [ <(> 1 {8)<h(8)erd8= f h{y)h{x-y)dy, 
2m J\ w \ =c Jo 

to evaluate the Laplace transform of the product 



n£i(«-«i) 
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Lemma 7. We have 

I 

xi...x n \ I dxj = 7 

'E?*i*j<i 7=i 



/■ 2fc T2k+n 

(3.20) A-a.^/'-^n^-p^- 



This lemma can be proved in a straight-forward manner by induction. 

4. Proofs 

We now give the proofs of our identities for the leading terms of the moments of the 
derivatives of A and Z. We begin with the proof of Theorem 2 for Z as it is slightly easier. 

Proof of Theorem 2. A differentiated form of the second formula of Corollary 1 gives us 



Uj- I \\Z A {e-^)Z A ,{e a ^)dA 

j=l 3 JU{N) h=1 



fc(fc-D L ^TT d e 2 ^ J= i J /■ , A 2 M 



--L 1 1 Kj<2fe J — 1 



+0 (jv* a+2fc - 1 ), 

provided that otj = otj(N) < 1/Af. Notice that 



and 



d 



(4-3) Jl^e")^ = = (-1)^(1)- 

So, 

(4.4) / \Z' A (l)\*dA N 



U(N) 

fc(fc+i) 



2fc , _ily-2fc 

d e 2 2-,=!^ /• a 2 w 



D^TT-- / e^=^ =^ T\dw 



rfojj k\(2iri) k J\ Wi \=i L -'-' (<''/ — n 



j=l l<j<2fc j 

rfc 2 +2fc-l\ 



,=i 



a=0 



+0(iV* 

The sign here arises as the (— l) kN from ()4.3|) cancels the same factor in (|4. lj) . we have a (— l) fc 
from (|4.2|) and we pick up the (— 1) 2 in (|4.1j) through writing the factor n^-y^i — w j) 
in ()3.7|) as A 2 (w) above. 
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To separate the integrals, we introduce extra parameters Li and move the Vandermonde 
polynomial outside the integral as a differential operator, getting 



(4.5) { _ x) ^f^^^/dL)e-T^ 



(4.6) 



+0(N k2+2k ~ 1 ). 
Next, we observe that 

d 



k\(2m) 



>Ej=i LiWi 



a=0,Li=N 



N 



da Ui<i<k( w i - °) 



a=0 



N 



so that ()4.5|) equals, without the O term, 



(4.7) 



. Hk+i) A 2 (d/dL) 

1) 2 i— 1 1 

; jfe!(27rf 



l) r ^j=i«j 



fc j JV 

ic. 2 



2 A: 



Li=N 



Introducing another auxiliary variable i, this can be expressed as 



(4 



d_\ 2k -Nt/2 



^X\d/dL)(i) M e 



3 Z)i=i LiWi+t/wi 



,2k 



km n i= i^ i= i 

This allows us to separate the integrals and we get 



Li=N,t=0 



k(k+ i) A 2 (d/dL)(d/dt) 2k e~ m / 2 k / 1 r 



(4.9) 

The integral evaluates to 
(4.10) 

as noted earlier. Thus, 



w 



2k 



-dw 



Li=N,t=0 



Li /2fc-i(2vZ^0 
(Lit) fc -V2 



(4.11) 



C/(JV) 



fc(fc+i) 

•1) 2 



^(i)r^= 

A(d/dL)(d/dt) 2k e- m / 2 f ^ L 2k - 1 I 2k . 1 (2^l~t) 



k\ 



n 

i=l 



(Lit) k ~y 2 



Li=N,t=0 



+0(N k2+2k - 1 ) 



So, letting 
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we have, by Lemma El that (|4.11j) equals 



2k 



kxk 



(4.13) 

Now we see, from (j'i.lOp . that 
(4-14) f t (L) = £ 

so that if fi < 2k — 1, then 



+ 0(JV' 



fc 2 +2fe-lx 



t=0 



r=0 



r!(2fc - 1 +r)! : 



(4.15) / M (^) = E 



r=0 



r\(2k- l-fi + r)\ \t 



L 



(2Jfc-l-/i)/2 



hk-i-fi(2v Lt). 



Therefore, (|4.13jl equals 



-l) a : , 

at 



2k 



e-^det , ( 

kxk \\t 



hk+i-i-AzVM)) 



■ 



t=0 



(4.16) +0(iV fc2+2fe - 1 ). 

Clearly detfc(aij) = det^afc+i-j^+i-j), therefore ()4.16|) can be written as 



(4.17) 



fc(fc+i) / d 
3) 2 { Jt 



2k x x AA (i+J-l)/2 

e^/ 2 det ( ( - 



t=0 



If we substitute x = iVt, then rf/rft = Nd/dx and we get 



(4.18) 



(-l)^^iV 2 



2k 



e~ x/2 det 

fexfc V \ x 



2\ (i+j-i)/a 



(-1) - S - JV 



dx 



2k 



J i+i _ 1 (2^) 
e- a ' /2 x- fc2/2 det (j i+ --_ 1 (2Vx) 



x=0 



x=0 



since detfc(M i+J 1 aj J ) = M fc2 detfc(ajj) as is seen by factoring M J out of the jth column and 
M i_1 out of the ith row. This completes the proof of Theorem 2. 

Proof of Theorem 1. Turning to Theorem l's proof, we begin as before, but with a 
differentiated form of the first formula of Corollary 1: 

2k k 



(4.19) f[jL / f\A A (e- ak )A A .(e a ^)dA N 



J=1 ««j h=1 

2fc 

fc(fc-i) 



d 



3= 

+ 0(N k2+2k - 1 ) 



\ daj k\(2m) k J\ w .\ = i 



A 2 (w) 



n 



i<i<fc (to, — a 



' i<i<2fc 



Y[dwj 



da 



d 



=° ds 



xAaWU = -a^(i) 
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provided that atj <C 1/iV. Now 

(4.20) 
and 
(4.21) 

hence setting a = 0, (j4.19|) becomes 
(4.22) 
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\a=0 



|A^(l)| 2fc ^ 



t/(A0 



2 A- 



^ 11 ^ A;!(2vrz) fc 



|Wi| = l 



FT kk* fio,- — a,-) 

LL l<i<2k K V 3=1 



+0(N k2+2k ~ 1 ). 

Introducing variables Lj as before, the above equals, without the O term 
k(k + l) Uf^d/da^A^d/dL) 



(4.23) (-I)" 



r2fc m/j. \^j/jn , pE^iLiWi-Nat) k 



k\(27ri) k 

Performing the differentiations with respect to the aj leads us to 
(4.24) 



Yldwj 



ctj=0,Li=N 



, k(k+i) A 2 (d/dL) 
-1) 2 



k\{2iii) 



Now we write 



IB; =1 



k / k 1 x fc 



3=1 



Li=AT 



Introducing the auxiliary variable t, (|4.24jl can be expressed as 
(4.26) 

^ ^ A\ Arfc _ h A 2 (rf/^)(|)^e-^ r e E<U * 



/i=0 



i .,, n /,r 



>../■ L. //A ^ l: _ h tf{d/dL)(d/dt) k+h e- Nt yr ( 1 

11 I 9tt 



fc(fc + l) ^-^ / fc \ _ r 



/i=0 



fc! 



8=1 



2tu 



\w\=l 



W 



21, 



-dw 



Li=N,t=0 
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Proceeding as before we arrive at 



(4.27) / \ti A (l)\ 2k dA N 

Ju(N) 



IU(N) 

rfc 2 +2fc-l\ 



x=0 



+0(N k 

Proof of Theorem 3. We now give the proof of Theorem 3. We rewrite equation (|4.1j) as 



2 /b j /i k 

(4.28) IT — / T\Z A (e-«»)Z A .dA N 

H d0l i Ju ^ h=l 



_i daij k\(2ni) k J\ w \ =1 H i<i<k (wi - ocj) ±Jr Wi 



+0(iV fc +2fc " 1 ). 

Introducing variables Lj as before, we can rewrite the main term above as 



(4.29) 



_ 1) M^il +feiV 2k d ^ / d \ « / d \ « ( i /■ e *,« rfu7 



n^-^W' n A n ^ 



fc! daj \dLJ f \ \dLiJ Al ^vrz 7 H=1 rg= x (w - or,-) ™ 

Now, by Lemma 6, the integral is 

(4.30) / TT & x 

Letting the variables in the ith integral be Xij we may express the product of the k integrals 

as 

(4.31) I ... I e E*=iE^i*««j Yl dxij. 

Ki<k 



l<j<2k 



We incorporate the factor e 2 ^j=i * into this product and have 
(4.32) / .../ e Zr=MzL^-m) TT dx 



Ej=l ^YlfLl x k,j<Lk l<i<k 

l<j<2k 
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We differentiate this product of integrals with respect to each ctj and set each ctj equal to 
yielding 

2k / k 



<« 3 > /„. ..-L IKE A ;1I 



ET=i x ^<Li J T,jii^k,j<L k j=i 



i=l ' l<i<fc 

l<j<2k 



We want to compute this integral by multiplying out the product and using Lemma 7. A 
good way to think about this is as follows. By equation ()2.18|) 

(4.34) {A x H V A k - A) 2k = ^ ( 2 £) {-A) m °A™ ■ ■ ■ A\ 

m£P* +1 (2fc) 

When we multiply out the product we will have a sum of [k + l) 2k terms, each term being a 
product of some number of factors {—N/2) and Xij. Let m £ PQ +1 {2k) represent a generic 
term in which {—N/2) appears m times, and factors xxj appear for mx values of j, and 
x 2 ,j for m 2 values of j and so on. When we apply Lemma 7 to this term, when we perform 
the integration over the variables xxj the answer is solely determined by mx, the number 
of different xxj that appear in this term. Therefore, we find that the product of integrals 
evaluates as 

2k\ / jV\ m ° ^2fc+mi j2k+m k 



{2k + mx)\ {2k + m k )\ 



s C)B) 

We now have that the quantity in equation (|4.29j) is equal to 
(4.36) 

( _ i)M ^_i) d k , d . _ , N . mo L 2k +mi L 2k + m k 



k\ ^ydLj^-ydLj 5! UJl 2 ) {2k + mx)\"' {2k + m k )V 

4=1 m<LP* +1 (2k) 

Now we need to carry out the differentiations with respect to the and set the Lj equal to 
N. We perform the differentiations \~[ k =1 d/dLi and obtain 

(4.37) 

(_1)=-t^+*jv 2 {d\ X ^ (2k\ ( N\" l ° Lf- i+mi L z k 



k\ \~dLj \^ \m) V 2/ (2fc - 1+mi)! "' (2fc + m fc )!' 



me?g +1 (2jt) 

Now the sum over nix, ■ ■ ■ m k is a symmetric function of the variables L{. Therefore, we can 
apply the second part of Lemma to obtain that the above, evaluated at Li = N is 



(4.38)/ \Z' A {l)\ 2k dA N 

JU(N) 



'U(N) 

/OU\ / AT\ m / AT2k+l+m i -i~_j \ 

^ fc 2 +2fc-l^ 



^ E )+0(jV , 



*;(fc+i) 

meP£ +1 (2fc) 
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which we rewrite as 
(4.39) 

m6P* +1 (2J;) 

Here the signs work out as in (|4.4jl . We factor l/(2k — i + rrii)\ out of the ith row. The 
remaining determinant has ith row 

(4.40) 

fc-i 

1, (2k-i + mi), (2k-i + m i )(2k-i-l + m i ), Yl^k - i - j + 1 + m<) 

i=i 

This determinant is a polynomial in the m, of degree + 1 + - • • + (A; — 1) = kik — l)/2 
which vanishes whenever — m, = j — i; moreover the part of it with degree /c(fc — l)/2 is 
precisely A(mi, . . . , m^) = Yli<i<j<k( m j ~ m «)- Consequently the determinant evaluates to 

(4.41) Yl (mj - mi- j + 

l<i<j<k 

This concludes the evaluation of b' k . 
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5. Numerical evaluation of b k and b' k 
We have the following values for bk- 

— 5 

bo = 



h = 



2 5 • 3 2 • 5 ■ 7 

3 ~ 2 7 • 3 4 • 5 2 • 7 2 • 1 1 

2275447 

2 18 - 3 10 • 5 4 • 7 3 • 11 • 13 

3700752773 
2 26. 314.56. 74. 112 . 132 -17-19 

3654712923689 
2 39 . 3 19 ■ 5 9 • 7 6 • ll 3 • 13 3 • 17 ■ 19 • 23 
53- 13008618017- 143537 



69 = 
ho = 





2 50 


. 3 28 . 5 13. 78. H5. 134. 172. 192. 

41 • 359 • 5505609492791 • 3637 


23 




2 68 


.335 


. 516 . 711 . 116 . 135 . 173 . 192 . 23 . 
757 • 45742439 • 60588179 • 13723 


29- 


31 


2 84 . 


3 42 


. 521 . 714 . us . 136 . 174 . 193 . 232 
652071900673 • 241845775551409 


•29 


•31 



612 = 



'13 



'11 



'15 



2 i05 . 355 . 525 . 717 . nio . 138 . 175 . 194 . 233 .29-37 

1318985497 • 578601141598041214011811 
2 i2i . 3 64 . 531 . 719 . 1112 . 139 . 177 . 196 . 234 . 29 2 • 31 2 • 37 • 41 • 43 

113 • 206489633386447920175141 • 51839 • 14831 
2 iso . 375 . 537 . 723 . 1115 . 1312 . 177 . 197 . 235 . 293 . 312 . 37 . 41 . 43 . 47 

4670754069404622871904068067089635254838677 
2 i74 . 390 . 542 . 72s . 1117 . 1314 . 1710 . 199 . 236 . 293 . 313 . 372 . 41 . 43 . 47 

107 • 194946046688455595346779341 • 996075171809335069 
2 203 . 3103 . 550 . 731 . 1120 . 1317 . 1712 . 1910 . 237 . 294 . 314 . 372 .41.43.47.53 

29547975377 • 3981541 • 1807995588661527603489333681461 • 1584311 
2 230 . 3117 . 557 . 737 . 1122 . 1319 . 1714 . 1912 . 239 . 295 . 315 . 373 . 4i 2 . 432 . 47 . 53 . 59 
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We have the following values for b' k : 

1 



K = 
b> 2 - 



2 2 -3 
1 



2 6 • 3 • 5 • 7 

3 2 12 • 3 2 ■ 5 2 ■ 7 2 ■ 11 

u 31 

4 2 20 - 3 10 • 5 4 • 7 2 • 11 • 13 

y 227 

5 2 30 -3 12 -5 6 -7 4 -ll-13 2 -17-19 

67 • 1999 

6 ~ 2 42 -3 19 -5 9 -7 6 -ll 3 -13 3 -17-19-23 



„ 43-46663 
b 7 - 



b' 



b[- 



256 . 328 . 5 13 . ? 8 . U 4 . 13 3 . 1? 2 . 19 2 . 23 

46743947 

2 72 . 334 . 516 . 711 . n e . 13 4 . 17 3 . 192 .23-29-31 

19583 ■ 16249 

9 2 90 • 3 42 • 5 21 • 7 14 • ll 8 • 13 6 • 17 3 • 19 3 • 23 2 -29-31 

, 3156627824489 

10 - 2 no . 3 55 • 5 25 • 7 17 • ll 10 • 13 s • 17 5 • 19 4 • 23 3 -29-31-37 

59-11332613-33391 
11 ~ 2 132 • 3 63 • 5 31 • 7 18 • ll 12 • 13 10 • 17 5 • 19 5 • 23 4 • 29 2 • 31 2 -37-41-43 

241 • 251799899121593 

0i 



12 2 156 - 3 75 - 5 37 - 7 23 -ll 15 -13 12 -17 8 -19 7 -23 4 -29 3 -31 2 -41-43-47 
, 285533 • 37408704134429 

13 ~ 2 182 • 3 90 • 5 42 • 7 28 • ll 17 • 13 14 • 17 10 • 19 s • 23 5 • 29 3 • 31 3 • 37 2 -41-43-47 

197 • 1462253323 • 6616773091 
14 ~ 2 210 • 3 100 • 5 50 • 7 31 • ll 20 • 13 17 • 17 12 • 19 10 • 23 7 • 29 4 • 31 4 • 37 2 • 41 • 43 • 47 • 53 

1625537582517468726519545837 
15 2 240 • 3 117 • 5 57 • 7 37 • ll 22 • 13 19 • 17 14 • 19 11 • 23 9 • 29 5 • 31 5 • 37 3 • 41 2 • 43 2 • 47 • 53 • 59 
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